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EQUATIONS FOR ISENTROPIC AND PLANE SHOCK 


FLOWS OF MIXTURES OF UNDISSOCIATED 

PLANETARY GASES 

By Victor L. Peterson 

Ames Research Center 
Moffett Field, Calif. 

SUMMARY 


Equations for aerodynamic quantities of interest in equilibrium isen- 
tropic one -dimensional flows and for properties across plane shock waves in 
arbitrary mixtures of planetary gases are derived under the assumption that 
the gases obey the perfect gas law. This assumption restricts the range of 
applicability of the equations to conditions for which the onset of chemical 
dissociation or ionization has not been reached although no requirement for 
the gases to be calorically perfect is imposed. Analytic functions used to 
specify the variation of specific heats with temperature are consistent with 
the assumption that molecular vibrations behave like harmonic oscillators and 
are independent of rotational degrees of freedom* Simplified equations giv- 
ing approximate results over more limited ranges of the variables are also 
presented. 

Sample comparisons of the aerodynamic behavior of a number of pure gases 
are made. It is shown that while some of the flow quantities characterizing 
aerodynamic flows are very sensitive to gas composition, the stagnation pres- 
sure coefficients for a number of gases differ by only a few percent from one 
gas to another. 


INTRODUCTION 


Considerable attention is now being directed toward the study of missions 
involving flight in the atmospheres of planets other than earth. The chemical 
compositions of the gases comprising these atmospheres are not known with any 
degree of precision at the present time. However, the presence or absence of 
certain constituents has been tentatively established for the near planets, 
and in some cases, abundances have been estimated. For example, evidence 
exists showing carbon dioxide to be present in the atmospheres of Mars and 
Venus in considerably larger amounts than found in the earth's atmosphere 
while methane and ammonia have been observed in the atmosphere of Jupiter. 
Various other gases such as nitrogen, argon, helium, and hydrogen have been 
considered as possible major constituents of planetary atmospheres. (See 



reference 1 for a survey of information on planetary atmospheres.) Thus, a 
number of gas mixtures considerably different from air are being associated 
with other planet atmospheres. 

In order to conduct analytical and experimental studies involving these 
proposed planetary gas mixtures it is necessary to understand the basic flows 
of these gases. Furthermore, if wind tunnels originally designed for opera- 
tion with air are converted to operate with various planetary gas mixtures, 
methods will be needed for relating aerodynamic performance in arbitrary gas 
mixtures to that in air. Many of the proposed planetary gas mixtures are 
unlike air in that they cannot be treated as perfect gases even at low temper- 
atures. Therefore, individual tables of flow information, such as those pre- 
sented in reference 2, have been computed for gas mixtures of interest. It is 
certainly not practical to prepare tables of aerodynamic parameters for each 
and every proposed atmospheric gas model so that the need for analytic methods 
is clearly evident. 

The present study is concerned with the development of equations for 
aerodynamic quantities of interest in equilibrium one -dimensional isentropic 
flows and adiabatic plane-shock flows for the ranges of conditions in which 
the gases considered obey the gas law p = pRT. The results are, therefore, 
generally applicable so long as the highest temperature in the flow does not 
exceed about 2500° K. This restricts the flow speed to the order of 2 to 4 
km/sec or below for most gas mixtures. The equations are derived for arbi- 
trary mixtures of monatomic gases such as helium and argon, linear molecule 
polyatomic gases such as nitrogen and carbon dioxide, and nonlinear molecule 
polyatomic gases not subject to internal rotation such as methane, ammonia, 
and water vapor. All of the polyatomic gases are admitted to be calorically 
imperfect with specific heats that vary with temperature. Within the frame- 
work of the assumptions results are obtained which are exact over the per- 
missible ranges of variables. Simplified approximate results applicable over 
more limited ranges of the variables are also presented. The present results 
that apply to real undissociated air may be compared to those in reference 3 . 
Unlike that study, however, specific heats are not restricted to vary in the 
first order only from their ideal gas values. 


SYMBOLS 


a 


A 

c p 

C v 

C ^stag 


sound speed 

area of flow cross section 

specific heat per mole at constant pressure 

specific heat per mole at constant volume 

pressure coefficient evaluated at stagnation conditions, 
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Naperian logarithm base 


energy per mole 
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enthalpy per mole 
molecular weight 
Mach number 

number of degrees of vibrational freedom 
pressure 

. . pV 2 

dynamic pressure, s ~— 

gas constant . — 

universal gas constant per mole 
number of atoms in molecule 
absolute temperature 

function of temperature defined by equation (8) 
function of temperature defined by equation (8) 
speed 

function of temperature defined by equation ( 15 ) 
mole fraction 

shock-wave angle measured from downstream flow direction 

C P 

ratio of specific heats , — 

angle of flow deflection across an oblique shock wave 
function defined by equation (22) 
molecular vibrational energy constant 
mass density 


3 


T ( ) ratio of ( ) for real gas to ( ) for ideal gas 

( ) partial derivative at constant volume 

(~ ) asymptotic value of quantity 

Subscripts 

i ideal-gas quantity 

j,k , Z indices referring to monatomic gases , polyatomic linear molecule 

gases, and polyatomic nonlinear molecule gases not subject to 
internal rotation, respectively 

n index denoting degree of vibrational freedom 

t total conditions for gas brought isentropically to rest (wind- 

tunnel reservoir conditions) 

air air quantity 


trans, rot, I translational, rotational, vibrational, and electronic energy 
vib, el J modes, respectively 

* quantity evaluated at sonic conditions 

1 conditions upstream of shock wave 

2 conditions downstream of shock wave 


ANALYSIS 


The analysis is divided into four sections: equations governing the 

thermodynamics of the gas mixtures are developed in the first, the second 
treats the one -dimensional isentropic flow situation, the third is devoted to 
plane shock flows, and the last presents simplified approximate forms of the 
equations. The equations presented only* apply to mixtures of thermally 
perfect undissociated gases. 


Thermodynamic Properties 


The thermodynamic properties of a mixture of gases can be calculated 
from the thermodynamic properties of the individual components or species. 
Thus, attention is first directed toward a general definition of the 
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necessary properties of an arbitrary pure gas- The ideas are then specialized 
to gases having specific types of molecular structure. Finally, equations for 
arbitrary mixtures of these gases are developed and a single equivalent gas is 
defined which has the exact thermodynamic properties of the actual gas mixture 
being considered. 

General considerations . - The total energy possessed by a gas molecule can 
be written generally as 


® “ -^trans + E int 


where E-t3 rans denotes the energy of translation of the molecule and Ep n t 
denotes the energy associated with the internal structure of the molecule. 
The internal energy Ej_nt can arise from three sources: 

1. E ro t rotation of the molecule about its mass center 

2 . E v -qj vibration of the atoms in the molecule relative to 

each other 

3. E e ^ excitation of the electrons within the molecule 

This study is restricted to gases for which these three phenomena can be 
assumed, with good approximation, to be independent. Therefore, 


E - E trans + E rot + E vib + E eZ 


The translational degrees of freedom are fully excited at all tempera- 
tures and quantum mechanical considerations show the rotational degrees of 
freedom to be fully excited at temperatures above a few degrees Kelvin. Thus, 
values for Etrans and E ro t are readily obtained by using the law of equi- 
partition of energy which states that each fully excited classical degree of 
freedom in a molecule contributes (l/2)RT to the specific energy. The energy 
associated with each degree of vibrational freedom in a molecule is determined 
from quantum mechanical theory by considering the vibrations to behave like 
harmonic oscillators. The energy in one degree of vibrational freedom is 

given by [RT(e/T) ]/(e©/ T - l) where 6 is the characteristic temperature of 
vibration. For a gas with K degrees of vibrational freedom 


■^vib " ^ 


I 


n-i 



1 


Note that each fully excited degree of vibrational freedom (T » G n ) behaves 
as two classical degrees of freedom by contributing RT to the specific 
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energy. The energy associated with electronic excitation is neglected in 
this study. At temperatures below those for molecular dissociation the 
electronic energy is small compared to the total energy. 

The specific heat at constant volume is defined as 


cv 



so that 

c v 1 /^ E trans\ 1 /dE ro A 1 fdE v ±i 0 \ 

X ~ R V dT ) v ) Y R V <W V 

It is apparent from the above equation that each degree of translational and 
rotational freedom contributes l/2 to the total c v /R, independent of temper- 
ature, while the vibrational degrees of freedom each contribute between 0 and 
1 to the total cy/R, depending upon the temperature. Thus, c v /R can be 
written in terms of a part independent of temperature (c v /R)^ and a part due 
to vibrations (cy/R)^*^ which is dependent on temperature so that 


c 


v 


R 




The relation between specific heats, Cp 


Cy = R, is used to obtain 



An ideal-gas value of the ratio of specific heats is defined as 


_ (c p /R) ± 
7 i = (cy/RJi 


(Cy/R)j. + 1 _ + 1 

(c v /R) i " + (c v /R) ± 


so that, finally, 
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1 + 


7, = 


(c v /R) . 




c v 

R 



C P 7 i . AbA 

R 7 . - 1 V R7 .. 

' 1 x ' vib 




( C p/ R ) 7 ± + (7 ± " l)(c v /R) vlb 
(c v /R) = 1 + (7 ± - 1 ) ( c v/R) vib 


(1) 


The enthalpy is defined as 


H = R 



T 



dT 


(2) 


This definition implies that the enthalpy is zero when the gas is in molec- 
ular form at a temperature of absolute zero. 

All of the thermodynamic properties needed in the analysis to follow 
have been written in terms of (c v /R)p and (cy/R)^^. Equations for each of 
these quantities will now be developed for each of three types of gases. 


Monatomic gases . - These gases have no molecular structure since they are 
comprised of single atoms. Thus, they have no energy in either rotational or 
vibrational degrees of freedom. Their entire energy is included in the three 
degrees of translational freedom so that 


= 3 

2 


V\ 

R/ . 

' vib 


- 0 


(3) 


y± 


= 1 + 


3/2 


5 

3 
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Polyatomic gases com prised of linear molecules . - The atoms comprising the 
molecules of these gases are arranged in a straight line. Each molecule 
possesses three degrees of translational freedom, two degrees of rotational 
freedom and one or more degrees of vibrational freedom depending upon the num- 
ber of atoms in the molecule. These facts lead to the following results. 



*\ 





= i + 


5/2 


7 

5 


( 4 ) 


Polyatomic gases comprised of nonlinear molecules not subject to 
internal rotation .- The molecules of these gases differ from the molecules of 
the linear molecule gases in that all atoms are not arranged in a straight 
line. Under these circumstances, energy can be stored in three, rather than 
two, degrees of rotational freedom. Furthermore, the number of vibrational 
degrees of freedom is given by n = 3s-6 instead of n - 3 S “5 as in the 
previous case . 

-\ 



J 


Mixtures of nonreacting gases . - The thermodynamic properties of a mix- 
ture of nonreacting gases are given by the general relations of equations 
(l) and (2) provided the functions (cv/R)^ and (cy/R)^^ are suitably 
defined. Let 
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Observe that (c Y /R) vib J^ = 0 and that ( 7 .^ = 5/3? (7 ± ) k = 7/5? and 

( 7 .) 7 = 4/3. Substituting these values into the equations and letting 
1 b 

(c v / R ) vi b = U ( T ) S^es 


where 




'V 


A k l 


~f) s U ( T ) = / x k u d T ) + / x Z u j( t ) 


Aa/T 




n=i 


«kn/ T 

e - 1 


3S-6 


u- 


;(T) = 


'In 


HJ ' 1 


n=i 




(6) 
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The enthalpy of the gas mixture, from equation (2), is 


H 




+ u(T) 


dT 


Let 


where 


H 



R 


T 

u(T) dT 



k 



3S -5 


V- 


k 


v — i 0, /T 

m T 


=o e 


n=i 


3S-6 


and v^T) = 7 


<V T 


n=i 


0, /T 

e - 1 


then 


H 

RT 


y a 

i-T + v(T) 

7.-1 

' 1 


(7) 


Note that the function v(T) is nothing more than Evik/RT. 
The gas constant R for the mixture is given by 


R = — 
m 


and the molecular weight m is given by 


m = ^ x jmj + ^ x k m k + ^ x z m z 
j k l 
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Summary of thermodynamic relations for gas mixtures. - 
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Figure 1.- Variation of the function u(T) 
with temperature for several gases. 



Figure 2.- Variation of the function v(T) 
with temperature for several gases. 


The functions u(T) and v(t) have 
been evaluated for a number of the 
gases considered in the compositions of 
planetary atmospheres. The results of 
these calculations are presented in 
figures 1 and 2. Values of the charac- 
teristic temperatures of vibration for 
a number of gases are given in table I. 
With the exception of air, the values 
given are consistent with those tabu- 
lated in references 4 and 5- The char- 
acteristic temperature used for air is 
that suggested in reference 3 . 


One -Dimensional Isentropic Flows 


The equations for quantities of 
interest in one -dimensional isentropic 
flows of nonreacting gas mixtures are 
developed in this section. The results 
of the preceding section showed that 
an equivalent gas having the thermo- 
dynamic properties of a specified mix- 
ture of gases could be defined. The 
equations of this section will be 
derived in terms of the thermodynamic 
functions of this equivalent gas. The 
thermodynamically equivalent gas has a 
ratio of specific heats 7 that is 
variable with temperature at high tem- 
peratures but which approaches a con- 
stant value at low temperatures. Flow 
properties for the gas mixtures com- 
puted using this' constant low tempera- 
ture value of specific heat ratio 7 . 
and the ideal gas equations such as 
found in reference 6 are called ideal- 
gas values. It has been found that 
equations for the ratio of quantities 
for the real gas to the same quanti- 
ties for the ideal gas take particu- 
larly simple forms. Thus, equations 
are presented for these ratios. Sim- 
plified forms of some of these equa- 
tions that give results accurate to 
within specified limits over limited 
ranges of the variables are derived in 
a later section. 
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The equations governing the four unknowns p, p, V } and T can he 
written as follows: 


Mass 

Entropy 

Energy 

State 


d(pVA) = 0 


dp c v dT n 

p " R T 


V dV + dH = 0 


p = pRT ; H = H(T) ; ^ 1 (T) 


( 9 ) 


The enthalpy and specific heat functions are given explicitly by equations (8). 

A relation between flow speed and temperature can be obtained directly 
by integrating the energy equation. The result is 


= 2RT 


?i 


1 7 . - 1 

/ 1 


2* - 0 + 71 ' 1 


y 


h v(T t ) - v(T x ) 


( 10 ) 


Both sides of equation (10) can be divided by the square of the sound speed 
to obtain the following relation for Mach number 


Mn 


V* 17 T. 


vy± - i) 


* - 1 1 + 7 - 1 — - 
y i 


L T 1 


v(Tt) - v(T,) 


( 11 ) 


It is recognized that by rearranging terms in equation (ll) and by adding 
unity to both sides the following result can be obtained 


T T 1 /T t ~ 


T k T t 

(Ti/Tt)i 


TiAt 



h 




v(T t ) 





( 12 ) 
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where use has been made of the fact that 





-l 


This function and all other ideal -gas functions employed in this analysis are 
given in reference 6. Equation (12) gives the correction to be applied to 
the ideal-gas results to account for the real-gas caloric imperfection. This 
correction is written in terms of Ti/l-t anc ^- either Tx or T-^ but it can be 
related to flow Mach number or speed in the following manner. In a flight 
application the static temperature Ti would be known and in a wind-tunnel 
application the total temperature T^ would be known so that the independent 
variable in equation (12) is the ratio Tj^/T-t* By use of the definition of 
t T /t 9 "k* 16 following equations for Mach number and speed can be obtained 






(13) 


2 2 7l RT 1 A Ti /T t \ 


(14) 


The right-hand sides of equations ( 13 ) and (l4) are functions only of T^T^ 
once either T x or Tt is specified, how that flow speed and Mach number 
have been related to temperature it is sufficient to determine each of the 
remaining isentropic flow quantities in terms of the static-to -total tempera- 
ture ratio Ti/lt and either static or total temperature. 


An equation for the real-gas correction to static-to-total density ratio 
is obtained by integrating the entropy equation. The result is 


1 



Now, recognizing that 

1 

7.-1 



i4 



and defining 


T pi/pt 


pi/ p t 

(pi/pt) ± 


allows the real -gas correction to the 


density ratio to be written 


X/pt 



u(T) 


dT 

T 


The integral in the equation above can be evaluated analytically. The result 
is 


where 


'T x 


u(T) ^ = w(Tj - w(T t ) 


w(T) =2^x k w k (T) +^x z w z (T) 
k l 


3S-5 


w k (T) 


n=x 


6 


n 0 kn/ T 

kn e 

T ^ta/ 1 






(15) 


W' 


3S-6 


,(T) = 


n=x 


(9 

e l n e in 


t / e ln /T 

e - 1 


7 ( e ^/T n 

In I e - 1 


J 


The functions w k (T) and w z (T) need be evaluated only once for each pure gas 
after which the function w(T) for any possible mixture of the pure gases can 
be formed. For convenience, the w(T) function has been evaluated for a 
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Figure 3*- Variation of the function w(T) 
with temperature for several gases. 


number of pure gases and the results 
are presented in figure 3. Finally, 
the real-gas correction to the static- 
to -total density ratio is written 


vpt - Cx/O 


\ i -1 - w(T t )] 

(16) 


The real-gas correction to the 
static-to-total pressure ratio is 
obtained in terms of the temperature 
and density corrections by using the 
equation of state. The result is 

_ p i /p t r -\r 'n 
Tp l /P t“ (Pj./PtJi ‘ V'V’V VVPlJ 

r Vi -1 [w(Ti) - w(T t )] 
= {j’EJT t ) e 



The real -gas correction to the dynamic -to -total pressure ratio is 
obtained as follows 


= q i /p t 
‘h./Pt (q 1 /p t ) i 


px v i g /p t = yx(Pi/Pt )M i 2 

(P 1 V 1 2 /P t ) i = 7 i (P 1 /P t ) i (M 1 2 ). 


T ^/P t = 7 ± T Pi/P t 


(18) 


An equation for the real-gas correction to the flow area ratio A/A* is 
developed from the mass continuity equation. Integrating the first of 
equations (9) gives 

P 1 ^iA = p^_V*A*. 


where the * denotes conditions evaluated at the sonic throat. Then, 
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but 


so that 


_L = p */ p t f* _1_ 
A* P x /p t a M i 


p */ p t / y » T * _1_ 

Pi/pt 7 r i T i 


/ft_\ _ (p^./pt)j A t */ t t)j _i_ 

w/i (p 1 /p t ) i y (T 1 /T t ) i Mi 


_ a/ a* 

t a/a* = (a/a*) ± 


Tp */ p t / r T^/T t 7 # 

T P 1 /p t 7 T T X /T t 7 i 


(19) 


Real-gas corrections to other flow quantities can be formed using 
combinations of those presented. 


Plane Shock Flows 


The equations derived in this section relate flow quantities on either 
side of plane shock waves. As in the preceding section, the equations are 
cast in terms of the thermodynamic functions of the equivalent single gas 
having the exact properties of any specified mixture of gases. 

The equations governing the four unknowns p 2 , p 2 , V 2 , and T 2 behind 
normal shock waves can be written as follows: 


Mass 


Momentum 


Energy 


State 


P-L V 1 = P 2 ^ 


P 1 + P 1 V 1 2 = P 2 + P 2 V 2 2 


Hi + | V x 2 = H 2 + | V 2 S 


p = pRT ; H = H(T) 


( 20 ) 


The enthalpy function is given explicitly by equations (8) 
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A number of approaches to the derivation of the shock relations are valid 
although in some the algebra becomes cumbersome. It will be shown that the 
density and pressure relations contain terms involving the static temperature 
behind the shock so a logical approach would be to first derive the relation 
for static temperature ratio across a shock. In the interest of simplicity, 
however, the density and pressure equations are derived first and then are 
combined to obtain the temperature equation. 

Manipulation of equations (20) results in the following quadratic 
equation for the density ratio p 1 /p 



This equation has the solution 


7 ± ( 1 + - J ( 7 i Mi 2 - 7i ) 2 + 2( 7 . 2 - l) 7 l Mx 2 ^ v(T a ) - v(TjJ 


7 x M i 2 (7i + 1) 


( 21 ) 


Let 


r] = (7 Mi 2 - 7 a ) 2 + 2(7 - 2 - 1)7 Mi 2 -=r- 


1 J-i 


T, "I 

v(T 2 ) - ^ v(T x ) 


( 22 ) 


and note that 


(l + 7- i)M! 2 


2 _ 


P i A 2 + ( 7± - 1)M X 2 


then 


_ Pp/Pi _ + (7j - D7 x M x 2 

rp 2 /P i (p 2 /p x ). 7.(1 + 7lMl 2 ) - J) 


( 23 ) 
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is obtained by first 


The equation for the static pressure ratio P 2 /p x 
noting that the momentum equation can be written 


2£=i + 

Pi 


7 Mi 2 



Substitution of equation (21) for p x /p into equation (24) gives 


(24) 


Now 


p 2 (i + r^i 2 ) + jq 
p i (r ± + i) 



27 i M 1 g - ( 7i - 1) 

(y± + i) 


so that 

Pa/P-L (1 + r^ 2 ) + ^/rf 

T P 2 / p x = (p 2 /p 1 ) i “ Zy^ - ( 7± - 1) 


(25) 


The real-gas correction to the static temperature ratio is obtained by 
using the thermal equation of state and the preceding equations for density 
and pressure. The result is 


T 2 /Ti 7 ± (l + 1 + (7 ± - l)(l + 7^ i £ )<Vn 


T 2 / t i (T 2 /T 


i' i 




i M x 2 - (7i 


1) 


(7, - l)M g + 2 


(26) 


An iterative procedure is required to solve equation (2 6) for the correction 
to the static temperature ratio since the function rj contains the unknown 
T 2 . A simple, rapidly converging, method requires writing T 2 in terms of 
the correction factor as follows: 


T 2 = 


T T 2 /Tx 



X 


Then, set T m equal to unity to obtain a first estimate of T 2 . This 

value of T 2 can be used to obtain a first estimate of t] which, in turn, 
allows the calculation of Tip • The new value of t t 2 /t I s ^ en used to 
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provide a second estimate of T 2 . The iteration is complete when the differ- 
ence between successive values of T Ta / Ti is within acceptable limits. It is 

shown in the section on approximations to the equations that the iteration can 
be avoided completely if results accurate to within several percent are 
acceptable . 

The real -gas correction to the ratio of flow speeds across a normal shock 
is developed from the mass, momentum, and state equations. 


V 2 

Vx 


7xMx 2 


1 + 7l Mi‘ 



(27) 


The ideal-gas equation for (Va/V-Jp is 


V 2 \ 2 + M i 2 (7 i “ 1) 

'Vi M 1 2 (7 i + 1) 


Combining equations (27) and (28) and substituting a previous relation for 

P 2 /Pi gives 


Y s /Y x 7 p(l + 7 -lV) " ^ 

V/V 1 = (Va/VxJi = 2 7l + 7xMi S (7 i - D 


(29) 


All of the remaining normal shock quantities can be developed from the 
results given above. The real-gas correction to the sound speed is obtained 
as follows: 


a 2 = / 7 a T g 
a i J 7 X Tx 


SO 



__ a s/ a l 

V/ax (a 2 /a 1 ) i 



(30) 
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The Mach number ratio is given by inspection 


Ms/Ml t v 2 /v! 


%/Mx - (Ma/Ml)i - T a2/ai = ( T p 3/Pi )( Taa/ai ) 


(31) 


The ratios of total conditions across normal shocks are given by the 
identities : 


r pt 2 /pt 1 


p t i / p ti (r^O 0 


P t; 


/ p O± 


r P 2 /P t 




Pt p/P 


' P t 2 /Pt 1 


2 


(VO, 


V P ti 


(32) 


The latter of equations (32) is obtained by making use of the fact that total 
temperature is conserved across the shock wave. 


The equations developed for normal shock waves are also applicable to 
the study of oblique shock flows. Values of the ratios T p /p > T p /p > 

t T 2 /T 7 an ^ T &2/ a i for -°hlique s ^ock waves are given by equations (23) , (25)^ 


(26), and (30), respectively, provided Mi sin p is used instead of Mi and 
the static temperature upstream of the oblique shock wave is the same as 
upstream of the normal shock wave. The flow turning angle & is related to 
the shock wave angle 3 by the following equation obtained from considera- 
tion of mass conservation 


tan (3 - 5) 
tan p 



2 + (y i - l)M 1 2 sin 2 3 
T p 2 /p i (7 i + l)^ 2 sin 2 p 


(33) 


The Mach number immediately downstream of an oblique shock wave can be 
determined from 
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( 34 ) 


M 2 Pi ai 1 (Ti/T 2 ) ± (p i /p 2 )i 

M 1 P 2 a 2 sm(p - s) (jp 2 /p^)(ja 2 /a^) sin (P - &) 


In the latter two equations the density, temperature, and sound-speed ratios 
are those corresponding to the upstream Mach number Mi and the shock wave 
angle (3. When M 2 is determined from equation (34), the remaining normal 
shock expressions, equations ( 32 ), can be applied to oblique shock flows. 


Approximate Results 


Simplifications to a number of the equations of the preceding sections 
are suggested by a cursory analysis. It will be shown for isentropic flows 
that substantial simplifications can be made, without great loss in accuracy, 
when the static temperature is relatively low (room temperature or less). 
Fortunately, this condition exists in a number of practical problems. For 
example, in most supersonic wind-tunnel applications the stagnation tempera- 
ture is so limited that the test section static temperature is nearly always 
of the order of room temperature or less. Also, in the flight application, 
the static temperature of the ambient atmosphere at altitudes sufficiently 
high for practical supersonic flight is generally low. It will also be shown 
that the equations for plane shock-wave flows can be simplified by imposing 
the additional condition that V 1 2 /RT 1 » 1 which implies the thermal energy 
of the gas to be much less than the kinetic energy of the flow. In addition 
to the simplifications, some rules to indicate the ranges of applicability of 
the approximate equations also will be derived. 

Isentropic flows . - Consider first equation (12) for the real-gas correc- 
tion to the static-to-total temperature ratio. When the static temperature 
is low then v(Tj_) -> 0. For these conditions equation (12) takes the form 


T Ti/T t 



1 + 


y± 



( 35 ) 


A significant result deduced from this equation is the fact that when the 
static temperature is low it ceases to be a variable in the problem. This is 
made evident by recalling that the value of y approaches the value of 7 ^ 
as the static temperature is lowered. Under these circumstances 7^/7 -*■ 1 

and the correction to the temperature ratio is a constant for a constant total 
temperature flow such as found in a wind tunnel. The value approached by 
t T /T asymptotically as T x -* 0 while T^ is held constant is given by 
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T T X /T t 


7 ± ~ 1 . 

= 1 + v(T t ) 


( 36 ) 


Substitution of equation (36) into equation ( 35 ) gives 

7a 


1 — 


r T 1 /T t ~ 7i r Tl /T t 


( 37 ) 


An estimate of the error introduced by using equation (37) for r T j ^ 

in place of the more complicated, but exact, equation (12) can now be made. 
Define 


f T = 
t T 


fT i/ T t ~ T r/ T t 


as being the fractional error. Substituting from equation (37) gives 


7 i 

f T = — - 1 

T T y 


Solving for 7^ in terms of the fractional error leads to 


1 1 + fn 


( 38 ) 


Equation (38) gives the value of 7 1 corresponding to a given fractional 

error in t T 1 /T-( ; ‘ This value of corresponds to a specific value of 

static temperature which can be determined from the equation relating 7 to 

temperature (eqs. (8)). Once static temperature is determined then equations 
(13) and (l 4 ) give the flow Mach number and flow speed, respectively, above 
which equation (37) can be used for calculating without incurring 

an error larger than f Tl ^. Results of sample calculations illustrating the 

real-gas temperature ratio correction for pure carbon dioxide are presented 
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in figure 4. The curves depicting exact results were obtained from equations 
( 12 ) and (l3)- The approximate results given by equations (37) and (13) are 
shown to agree closely with the exact curves for Mach numbers where the static 



Figure 4-- Comparisons of exact and approximate results for the real gas correction to the 
static-to -total temperature ratio in an isentropic flow of carbon dioxide. 


temperature is low. The asymptotic result given simply by equation (36) is 
shown to be essentially exact. The predicted and the exact lines of constant 
percent deviation from the asymptotic values of are also shown. Note 

that the theory is conservative in predicting the Ma,ch number for which 
approximate results deviate a specified percentage from the exact results. 

The same general procedure used to estimate the range of applicability of the 
simplified equation for the real-gas correction to the static-to-total tem- 
perature ratio can be used to provide similar estimates for the corrections 
to the other isentropic flow parameters. 


In the case of the real-gas correction to the static-to-total density 

ratio, it is observed that the function e W ^ Tl ^ in equation (16) approaches 
unity as the static temperature T^ becomes low so that it is convenient to 
define 


l/p t ~ VTi/T t 




-w(T t ) 


(39) 
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Equation ( 39 ) gives the value approached by t / asymptotically as T 

p t 

becomes progressively lower while the total temperature is held constant. 
Using equation (39); equation (l 6 ) can be rewritten as follows 




r T 1 /T, 


XX 


t Vi' 1 - 


X/X 


Xx 


w(Ti) 


(40) 


An estimate of the flow Mach number for which a given fractional error is 
incurred by using the simpler equation, equation ( 39 ); in place of the exact 
result of equation (40) can now be made. The approximation for r f 


T 1 /T t 

given by equation ( 37 ) is substituted into equation (40) and the term e 
is dropped so that 


w 


(Ti) 


Now, write 


p xX 




(41) 


Expanding the right side in a Taylor series and dropping higher order terms 
gives 


l 



Define the fractional error as 


f r = 

T p T 


V p t - XX 


PxX 


and substitute equations (4l) and (42) into this definition to get 
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Solving for y x gives 


7, = 


i 1 + (7 ± - l)f Tp 


(43) 


Equation (43) gives the value of specific heat ratio corresponding to a given 
fractional error in t / resulting from using the simplified result of 

p i/ p t 

equation (39). By methods discussed previously, the value of y given by 

equation (43) can be related to static temperature and, in turn, to either 
Mach number or flow speed. 


An approach similar to that used above gives the following ratio of real- 
gas correction to the static-to-total pressure that is approached as the 
static temperature T 1 is diminished and the total temperature held constant 



t 



(44) 


so that 


F Pl/P t \ Tr 


2 j_ 

T i/ T t\ 7:L w(T x )- 


'Tx/Tt j 


Pi/p t 


The fractional error defined by 


f-r.. a 


r P x /P t - T P x /Pi 
T Px/Pt 


can be written 


fr P ^7 X " V 7 , - 1 


(45) 


so that the value of the specific heat ratio corresponding to a given frac- 
tional error in Tp /p^ resulting from using the simplified result of 

equation (44) is given by 


7 


l 



(46) 
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Plane shock flows . - Two levels of simplification to the plane shock 
relations are considered. The first of these eliminates the need for the 
iterative solution across the shock wave. The second involves obtaining the 
equations for the asymptotic values of the real-gas corrections approached as 
the static temperature is reduced while the total temperature is held constant. 

The equations for the plane shock flows are written in terms of the func- 
tion t) defined by equation (22). Evaluation of this function requires a 
knowledge of the static temperature downstream of the shock wave. This tem- 
perature is not known a priori so that an iterative solution of equation (26) 
for t t 2 / Ti is required in order to establish a value for T 2 and subse- 
quently tj. This difficulty can be circumvented, with only small loss in 
accuracy, by recognizing that T 2 is never more than a few percent different 
from the total temperature and that the r\ function is not very sensitive to 
errors in T 2 of this order. Because of this, T 2 can be replaced by T^ 
in the equation for tj . A comparison of exact results for the real-gas cor- 
rection to the static temperature ratio across a normal shock wave in pure 
C0 2 to those obtained using the simplified r\ function is made in figure 5* 
The differences between the exact and the approximate results are shown to be 
small. 



Mach number, M 

Figure 5*- Comparisons of exact and approximate results for the real gas correction to the 
static temperature ratio across a normal shock in carbon dioxide. 


The results of figure 5 also show that the real-gas correction to the 
static temperature ratio approaches an asymptote as Mach number is increased. 
This behavior is typical of all of the real-gas corrections to the shock 


parameters as it is for the corrections to the isentropic flow quantities. 
Equations for asymptotic values of all the functions can be readily obtained 
after finding the form of the t] function for T x -> 0 and T 2 = Under 

these circumstances, y ± -* 7 ^ so that 


T) = J ± 2 Mx 
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(r ± 2 - 1 ) 

T Ti/T t 


Substituting this relation for if into equations ( 23 ), (25), ( 26 ), and ( 29 ) 
and by dropping terms of order 1/7 Mi 2 = RTx/Vi 2 gives the following 
equations for the asymptotic values of density, pressure, temperature, and 
speed ratios, respectively. 



y± 


(7 i - 1) 

k - ^ 

J T T!/T t 


( 47 ) 


T 


P 2 /Px 




(7i 2 - 1 ) 
T Ti/T t 


( 48 ) 


T T 2 /Ti 


1 

2 


( y ± + 1) 

T Tx/Tt 



(7j 2 - 1 ) 
T Ti/T t 


( 49 ) 


T V 2 /Vx 





(7 ± g - 1) 
T Ti/T t 


( 50 ) 


The asymptotic value of the real-gas correction to the sound speed ratio is 
obtained from equation ( 30 ). The result is 


r a 2 / a x 



(51) 
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Equations for the asymptotic values of the corrections 
parameters are given by combinations of equations (47) 


to the remaining 
through ( 51 ). 


DISCUSSION 


Using the equations developed in this study it is a simple matter to 
determine the effects of changing gas composition on quantities of interest 
in isentropic and plane shock flows . Some sample comparisons of the behav- 
ior of a number of gases will be made and discussed. 

Consider first the situation of a body flying in an atmosphere. The 
variation of total temperature with flight Mach number and flight speed for a 
number of pure gases is shown in figure 6. These results were obtained by 
means of equations (12) , ( 13 ), and (l4) and are for an arbitrarily chosen 
ambient temperature of 200° K. It is shown in figure 6(a) that changes in 
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(a) Variation of total temperature with Mach number. 


Figure 6.- Effect of gas composition on the variation of total temperature with Mach number 

and speed; - 200° K. 
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gas composition have large effects on total temperature when results are 
compared at the same Mach number. However, as shown in figure 6(b), the com- 
parisons of results for the various gases are considerably different when made 



Speed, V| , km/sec 

(b) Variation of total temperature with speed. 
Figure 6.- Concluded. 


at the same speed. For example, the difference between the total temperature 
in carbon dioxide and that in nitrogen at comparable Mach numbers grows with 
increasing Mach number to exceed 800° K while the difference never exceeds 
about 150° K when evaluated at comparable speeds. In contrast to this result 
the differences between the curves for helium and air are larger when compared 
at the same speed than when compared at the same Mach number. In addition to 
magnitudes, the algebraic sign of the difference is not always the same for 
the results presented as a function of Mach number as for the results pre- 
sented as a function of speed. Considerations of this sort indicate that when 
comparing aerodynamic quantities for various gas compositions it will be nec- 
essary to exercise care in drawing conclusions since the conclusions can be 
strongly dependent upon the basis of the comparisons. 

As a second illustration of the aerodynamic behavior of various gases 
predictable by the equations derived in this analysis, consider an application 
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Figure 7.- Effect of gas composition on Mach 
number obtainable in a wind tunnel. 


Figure 8 .- Effect of gas composition on speed 
obtainable in a wind tunnel. 


to wind-tunnel flows. In many wind tunnels the total or reservoir temperature 
is held constant during a test as opposed to the flight case where the ambient 
temperature is fixed. For a total temperature of 2500° the ratio of flow 
Mach number obtainable in each of a number of gases to Mach number obtainable 
in air is presented in figure 7 as a function of air Mach number. Results 
presented in this form readily show estimates of the performance capabilities 
of presently operating air facilities when operated with planetary gas mix- 
tures. It is apparent from the results of figure 7 that the flow Mach number 
obtainable in a given nozzle is strongly dependent on gas composition. A noz- 
zle designed to produce a given Mach number in air will produce higher Mach 
number flows of gases having a smaller ratio of internal energy to total 
energy than air and lower Mach number flows of gases having a larger ratio of 
internal energy to total energy than air for a given reservoir temperature. 

The ratio of wind-tunnel flow speeds obtainable with each of a number of 
gases to that obtainable with air is shown in figure 8 as a function of air 
Mach number for a total temperature of 2500° K. The relationships between the 
various gases in terms of flow speed obtainable in a wind tunnel designed for 
a given Mach number in air are quite different from those shown in figure 7 in 
terms of Mach number. All the differences between the comparisons illustrated 
in figure 7 and those shown in figure 8 are due to differences in sound speed 
in the various gases. 
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As a final comparison of the aerodynamic behavior of a number of possible 
constituents of planetary atmospheres , consider the stagnation pressure coef- 
ficient. This quantity is of particular interest since a number of theories 
for predicting pressures and forces on entry bodies relate the predictions to 
the stagnation point pressure. The stagnation point pressure coefficients for 
a number of gases normalized to that for real air are presented as a function 
of flight speed in figure 9- It is interesting to note that while many of the 



Figure Effect of gas composition on stagnation point pressure coefficient; T 1 = 200° K. 


quantities characterizing aerodynamic flows are very sensitive to gas composi- 
tion, the stagnation pressure coefficients for a number of gases differ by 
only a few percent from one gas to another. 


CONCLUDING REMARKS 


Equations for aerodynamic quantities of interest in isentropic one- 
dimensional flows and for properties across plane shock waves in arbitrary 
mixtures of undissociated planetary gases have been derived. It was shown 
that the equations could be cast in forms not prohibitively complicated by 
making the following principal assumptions for the gas model: (a) all gas 

components of a mixture are thermally perfect so as to be governed by the 
perfect gas law p = pRT; (b) molecular vibrations behave like harmonic 
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oscillators with no coupling between rotational and vibrational degrees of 
freedom; and (c) electronic energy is negligible compared to the total energy. 
Simplified equations giving approximate results over more limited ranges of 
the variables were also obtained along with estimates of the error introduced 
by the approximations. 

Various uses for the equations were demonstrated by making sample compar- 
isons of the aerodynamic behavior of a number of gases. In one comparison, 
estimates of the flow Mach number and speed capabilities of a wind tunnel 
operated with a number of pure gases in terms of its performance in air were 
presented. It was shown for the same reservoir temperature that a nozzle 
designed to produce a given Mach number in air will produce higher Mach number 
flows of gases having a smaller ratio of internal energy to total energy, at a 
given temperature, than air and lower Mach number flows of gases having a 
larger ratio of internal energy to total energy, at a given temperature, than 
air. It is also shown that while some of the quantities characterizing aero- 
dynamic flows are very sensitive to gas composition, the stagnation pressure 
coefficients for a number of gases differ by only a few percent from one gas 
to another. 


Ames Research Center 

Rational Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 13, 1964 
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TABLE I.- CONSTANTS FOR SELECTED GASES 


(a) Linear molecule gases 


Gas 

Molecular 

weight 

Vibrational energy 
constants, 0 n , °K 

U) 

CO 

28 

3121.85 

C 0 2 

44 

959 - 9 ( 2 ), 1928, 3379-5 

h 2 

2 

6322.5 

NO 

30 

2738.44 

n 2 

28 

3394.3 

n 2 o 

44 

845 - 8 ( 2 ), 1853-1, 3218 

O2 

32 

2273-35 

air 

29 

3055-5 


(b) Nonlinear molecule gases 


Gas 

Mole cular 
weight 

Vibrational energy 
constants, © n , °K 

(1) 

ch 4 

16 

1879(3), 2195(2), 4192.1, 4344.7(3) 

h 2 o 

18 

2294.4, 5253, 5402.7 

nh 3 

17 

1340.08, 1392.58, 2341.6, 
4798.7, 4801 , 4911 


■‘‘Number in parentheses denotes the number of degrees of 
freedom to which the constant it follows applies. 
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